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Abstract 

We consider a birth-death process with the birth rates iX and death 
rates i/i+i(i—l)9, where i is the current state of the process. A positive 
competition rate 9 is assumed to be small. In the supercritical case 
when A > fi this process can be viewed as a demographic model for a 
population with a high carrying capacity around . 

The article reports in a self-contained manner on the asymptotic 
properties of the time to extinction for this logistic branching process 
as 9 — > 0. All three reproduction regimes A > fi, X < /x, and A = \i are 
studied. 
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1 Introduction 

One of the basic population models with continuous time is the linear birth- 
death process (Xo(t),t > 0) with fixed birth and death rates A and /i per 
individual. This is a simple example of a branching process describing a 
population of independently reproducing individuals having three different 
reproductive regimes: supercritical (A > u), critical (A = /i), and subcritical 

(a<a*). 
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The properties of the linear birth-death process Xq(-) and its time to 
extinction tq are well-known, see for example [3 pp. 270-2]. In particular, 



E m X (t) = me 



(X-fi)t 



and 



(r < t) = { V A-MeC-*)* 



, in the sub- and supercritical cases, 
(™) m , in the critical case, 



where P m (-) and E m (-) stand for the conditional probability and expectation 
given that the corresponding birth-death process starts from the state m. 
It follows that in the supercritical and critical cases E m ( T o) = oo and in the 
subcritical case E m (ro) < oo. Letting t — > oo one obtains the extinction 
probabilities 



t (r < oo) 



1, in the subcritical and critical cases, 
(^) m , in the supercritical case. 



Moreover, it is easy to see that in the subcritical case 
lnm + ln(l — X/ fi) + x\ ^ -(e~ x ) 



t < )->e (e >, m-^oo, (1) 

fi-X ) 

and in the critical case 

E m (to < mx ) — > ex P{ — i^ x )~ 1 }^ m^-oo. (2) 

The absence of competition among individuals is a major weakness of 
the linear birth-death population model. A natural modification of this 
simple-minded model is to introduce extra deaths due to competition. We 
consider an indexed birth-death process (Xg(t),t > 0) taking non- negative 
integer values i G {0, 1, 2, ...} and having time homogenous jump rates 

Fi(X e (t) = i + 1) = Xit + o(t), with Xi = iX, 

Fi(X e (t) = i-l) = mt + olt), with iii = in + i{i- 1)0, (3) 

Fi(X (t) =i) = l-(Xi + m)t + o(t), 

as t — > 0. The key parameters (A, fj,, 6) of the model are the birth, death, 
and competition rates providing the following description of the demographic 
dynamics until the process hits the absorption state i = 0. 

Given the current population size i > 1, the next change in the popula- 
tion size is caused either by a birth or by a death of a particle. It is assumed 
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that coexisting particles give birth independently of each other at rate A per 
particle, so that interaction among particles does not influence birth events. 
Particle death is modeled by two parameters: parameter \x gives the death 
rate per particle "due to natural causes" and parameter 8, usually assumed 
to be small, quantifies the death rate due to competition pressure (factor 
i(i — 1) appearing in front of 8 represents the number of pairs of competing 
particles). Putting 8 = brings us back to the linear birth-death process 
Xq(-) mentioned in the Introduction. 

The process Xg(-) is an example of the so called logistic branching process 
studied in [10] along with its continuous state counterpart. The birth-death 
framework allows for a more detailed analysis in this special case. The most 
conspicuous new feature of Xg(-) compared to the linear birth-death process 
Xq(-) is the existence of a threshold value 



in the supercritical case. Obtained from the equation Aj ~ Hi the threshold 
value ig splits the state space in two parts. For i < ig the process Xg(-) tends 
to grow while for i > ig it tends to decrease. A relevant biological interpre- 
tation of this threshold value is the carrying capacity of the environment for 
the population in question. 

In Section[2]we summarise some useful properties of the time- homogeneous 
birth-death processes. It follows, in particular, that the quadratic form of 
the death rate compared to the linear birth rate ensures that our birth-death 
process with competition goes extinct with probability one (in contrast with 
a supercritical linear birth-death process which never dies out with a posi- 
tive probability). One of the most interesting characteristics of the process 
Xg(-) is the random time to extinction Tg. 

If 8 is small, the competition component i(i — l)8 is much smaller than i[i 
for i <C 8 , so that the process Xg(-) at relatively low levels can be approx- 
imated by the linear birth-death process Xq(-) with parameters (A,//) and 
the same initial state Alo(0) = m. This is done using a coupling construction 
presented in Section [3j 

Section [4] presents the main asymptotic results for expected value and 
distribution of the time to extinction Tg as 8 — > 0. The remaining sections 
contain the proofs. 




A — fj, 



(4) 
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2 General properties of time homogeneous birth- 
death processes 

Next we give a short summary of useful results for a time homogeneous 
birth-death process with birth rates A, and death rates Hi-> some of these 
properties can be found in [6] and [7]. An important probability 



satisfies a recursion 



implying 



Using notation 



Qi = Pj (reach i + 1 before 0) 



Qi — , ~i~ , Qi — lQi 

Aj + Hi Ai + /Xj 



1 A, 1 



1 — Qi Hi 1 — 



7T = 1, TTj = f-, IIo = 0, life = > 7T 

Al • • • Aj 



J=0 



we derive Qi = rr 1- - 

U-z + l 

More generally, for i £ (k, n) C (0, oo) 
Pj (reach n before k) 
Pi (reach k before n) 



n—n fc 
n n - n fc ' 
n n -Ui 



n„-n fc 

Using this we can compute the conditional jumping probabilities 

Pi (first jump goes down | reach i + 1 before k) 

Hi Pj_i (reach i + 1 before k) Hi IIj_i — IL^ 



Aj + Hi Pj (reach i + 1 before fe) Aj + /ij IL; — II^ 
which in turn lead to the recursion 

Pi = Ej(time to reach i + 1| reach i + 1 before k) 



1 ( Hi H' i 'U-,, , ,/,-. 



A l + /ij Aj + ^ rij - 
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resulting in a difference equation 

' A,ill,, : II;/ ' 
which is easily solved as 



,. _ Uj - u k + ni(Ui-i-u k ) k p k = 1 



ok >}i W ~ ^k) 

Pi (n,-n fc )(n m -n fc ) ' u 



Similarly, for G the conditional jumping probabilities 

Pfc (first jump goes up | reach k — 1 before i + 

Afc Pfc+i (reach — 1 before i + 1) _ Afc Ili+i - U k +i 
Xk + fJ'k Pfc (reach A; - 1 before i + 1) \k + Hk n i+ i — 

give the recursion 

= Efc( time to reach A; — 1| reach fe — 1 before i + 1) 

= 1 + A fc iii+i - n fc+ i i+1 + 

resulting in a difference equation 

pi+i = n »+i ~ n fc + A fc n m - n fc+1 i+1 i+1 = i 

which implies 



J—K 



Observe that for all v > u > relations ([5]) and (|6| bring 

Eo« _ _ (n^+i - nj)(iij - n n ) 

Pi - 1^ Ph - A^-m^+i-nj ' { ' 

i=u+l k=u+l j=u+l 3 3y v+l u > 

This is a confirmation (in terms of the first moments) of the statement in 
claiming that the corresponding conditional hitting times are equal in 
distribution. 
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The expected absorption time is given by the formula 

i— 1 oo _. 00 n 

Ej(time to reach 0) = V vr,- V = V — — . (8) 

J \ k7Tk \ k7Tk 

j=0 k=j+l K K k=l 

Indeed, if we denote the last expectation by ai, then the following recursion 

1 Aj fii 
oti = t — ■ h t — ; Oi+i + t — , oti-i, 

Aj + Hi Ai + Hi A i + Mi 

takes place with ao = 0. From this recursion it is straightforward to derive 
formula It follows from Q that 

X j 

Ej(time to reach i — 1) = 7r^_i > . (9) 

In particular, for the subcritical linear birth-death process formula ^ 
gives 



k 

k=m+l 



^ \fe=l fc=m-t 

1 / m 00 

— { jV 1 + ln(l " *) + a - " £ 

\fc=l k=m+l 



S k 



with s = X/ /J,, implying 

„ , . Inm + 7 + ln(l — A/it) , s , 

E m r = 1 ; ^ + o 1 , m -> 00, 10 

M — A 

where 7 = 0,577... is Euler's constant. This complements the weak conver- 
gence ([!]) in terms of asymptotic equality of the corresponding expectations. 



3 A coupling to the linear birth-death process 

To partially extrapolate the nice properties of the linear birth-death process 
Xq(-) to the process with interaction Xg(-) one can use the following coupling 
construction (cf [1]). 

Consider a bivariate Markov process (Xg(-), Xq(-)) with transition rates 
given in the next list. 
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Type of transition (0 < i < j) Transition rate 

-)■(? + + 1) Xi 

(i, i) — > (i — 1, i — 1) iii 
(i,i) ->•(£ — l,i) - 1) 

-> (i+ Xi 

-»• (*,i + l) Aj 
->•(*, J- 1) MJ 

The process is constructed in such a way that -X#(i) < Xo(£) for all t > 0, 
and the marginal distributions of jfo(-)) coincide with those of -Xfl(-) 

and -Xo(-), respectively. 

An important question here is how long this bivariate process stays at 
the diagonal if (Xg(0), Xo(0)) = (m,m). Let Kg be the number of jumps of 
the process (Xg(-), Xq(-)) until separation, if the components stay together 
until extinction we put Kg = oo. We show below that 

™ / v (m + n)n9 

P(m,m)(«* < «) < S— ~ , (11) 

where P( mjm )(-) stands for the probability conditioned on the bivariate pro- 
cess starting from the state (m,m). 

Suppose 9 — > and the starting level m is fixed. In the subcritical 
and critical cases the total number of births and deaths in the linear birth- 



death processes is almost surely finite and due to (11) we may conclude that 
T e r o almost surely. Moreover, since a supercritical branching process 
conditioned on extinction behaves like a subcritical branching process, we 
obtain that tq — > tq almost surely provided tq < oo. This observation is 



summarised in the next section as a part of Theorem 4.1 



To prove (11) observe that 

Kg = inf{k : U k / V k }, 

where Vk)k>o ar e the consecutive states visited by of the process (Xg(-), Xq(-)). 
Note that the only way for the bivariate process to get off the diagonal is the 
move (i, «)—>•(£— 1, i) having the probability j^^g^fj which is negligible, 
if the current level i is not too high. Since 

m+n-l q,^ _ 

P(m,m)(«0 = n \ K e > n - 1) = ^2 — — _ F^ m) (U n -i = i) 

i=l 

< 9 F TT < {m + ^ 
A + n K ' X + \i 
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we derive ( 11 ) 



\m,m)( K 6 < n) = 1 - Y[ ^(m,m){Ke > k\Kg > k - 1) 



v (m + A;)0 (m + n)n9 

4 Main Results 

We claim that as 9 — > the following two limit theorems hold for the birth- 
death process defined by ||. 

Theorem 4.1 If Xg(0) = m, where m is a fixed positive integer, then 
(i) in the subcritical and critical cases when A < fi 

W(m,m)( T 6 -> t ) = 1, 
(wj m the supercritical case when A > fi 

^(m,m)(Te ->■ To|t < Oo) = 1, 

and for any x > 

P(m,m)fa > xciv^e C2/e )|r = oo) -> e -1 ', 

where 



ci 



A(A - VWm, c 2 = A-yu-/iln(A/At). (12) 



Theorem 4.1 (i) and the first part of (ii) are proven in the previous 



section. The proof of the second part of Theorem 4.1 (ii) is given after the 
proof of the first part of the next theorem. 

Theorem 4.2 If Xg(0) = mo and 6mg — > a > 0, then 
(i) in the supercritical case when A > fi 

E mg (T e )~ Cl Vee C2 / e 



with positive constants c\, C2 given by (12), and for any x > 
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(ii) in the subcritical case when X < fi 



Em fl (re) 
and for any x > 



Infae- 1 ) + In ^ + In 



n—\+ 



a+7 



fj, — A 



T(9 < 



\nia8~ 1 ) +ln^ + In 



s (-e— ) 



(13) 



/j — A 

('raj in £/ie critical case when A = /j 



(14) 



(tt/2) 3 / 2 



0jU 



The asymptotic formulae for K mg (Tg) in Theorem 4.2 (ii), (hi) are ver- 
ified by simulations as shown in Figures [TJ [2] Comparing the asymptotic 
formula (10) for the linear birth-death process to the that for the process 

a the average 



with competition (j 1 3|) we see that as 9 — > and 9mg 

[i — A + a 



survival time reduces by 

E me fa) -E me (T e ) 



1 



In 



[i — A n — A 



As one would expect, this difference becomes small for larger values of /x 
and/or smaller values of a. 

Notice that Theorem |4.2| (ii) provides with a counterpart of the weak 
convergence ([!]) for the linear birth-death processes, however, we could not 
find a counterpart of (|2j) in the critical case. The following lemma plays a 
crucial role in the asymptotic analysis of all three cases. 



rj-l n+6i 



Lemma 4.3 For our particular model the function ttj = Y\ i= o 
isfies the approximation 



sat- 



ttj = (l+jO/ri-We-wW/'il + ritf)), j > 1, (15) 

where W(x) = x — x In — [i In and for any fixed T > 

sup \r)j(0)\->0, 0-»O. 
i<j<T/e 
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Figure 1: Averages of 100 simulations for each value of the death rate are 
plotted against the values predicted by Theorem 4.2 (ii). Choice of parame- 
ters: initial population size mg = 1000, competition strength 6 = 0.001, and 
birth rate A = 1. 
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Figure 2: Averages of 100 simulations for mg = 1000, 9 = 0.001 are plotted 



against the values predicted by Theorem 4.2 (iii). 



11 



Proof Observe that 



3-1 



' 111 TTj = 6>E ln 

fc=0 

In 



o 



fi + kd 
A 

A 



fc = 



3-1 



k=0 r 



It remains to verify that 



sup 

i>o 



n 

fc=0 



fi + jd 



□ 



5 Proofs for the supercritical case 



In this section we first prove Theorem |4.2| (i) and then Theorem 4.1 (ii) 
borrowing key ideas from [I]. 

We start by considering the supercritical Xg(-) with the initial state 
Xq(0) = ig given by Q. It will take a geometric number Kg ~ Exp(l — Qi g ) 
of returns to the initial state from above before the extinction event. Let 
Tg be the time needed for Xg(-) to enter the level ig from above, and To be 
the absorption time counted from the last entrance moment to the state ig 
from above. If there were no visits of ig from above, we put r' e = rg and 
t'q = 0. Clearly, Tg is the sum of Tq, r^', and of Kg independent durations of 
the corresponding excursions. It follows that the statement (ii) of Theorem 



4.2 is a straightforward consequence of the next three lemmata. 



Lemma 5.1 In the supercritical case as 8 — > 



l-Q 



in 



A 3 / 2 



(16) 



Lemma 5.2 In the supercritical case the expected duration Mg of an ex- 
cursion starting from ig and returning to ig from above satisfies 



Me 



2n8 1 
A A — fi 
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Lemma 5.3 Under the assumptions of Theorem (i) for any fixed 
positive e 

E mg (rl ) ) = o(e^), E me (4') = o(e^ e ), 0^0. 



Proof of Lemma 



5.1 



According to Lemma |4.3| 



It is shown in Section Jij that 1 — Qi e = ■ 



7Ti. 



(17) 



and in view of IL 



X—fi 



we arrive at (16). 



Proof of Lemma 15 . 2L We show first that 



OO j 

~W 



1 



:V2^9e C2/e 



□ 



(18) 



by dividing the left sum in two parts using the threshold 
Observe that by definitions of i$ and 1$ we have Hj > 7r, e • 2 J "~ ' 9 for j > 
Thus applying (17) we obtain 



2A-u- 



1 

£ — = O(0/vr i9 ) = O(^/ e ). 



On the other hand, due to (15) 



j=l J 3 







fly/// 



w ^' e ds. 



As the function W(-) has its maximum over the integration interval at the 
point A — fi we conclude using the Laplace method that 



V2^9e c ^ e 



{ x i n 3 ( A - 



6 1 

~\7n 



completing the proof of (|18[ ). 
Combining ^ and (|9| we get 

M e = ^if, ^2 y~' where ^(i) = min 



k=l 



iLn m 
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Relations (fT7|) and (|18[) give 

/2tt6> 1 



, AfcTTfc V A A - fx 

and it remains only to observe that tpki^e) ~^ 1 uniformly over k larger than 
e/9 however small is a fixed positive e. 

□ 

Proof of Lemma |5.3[ Notice that due to Q 



>0 



E(tq) < Ej e (time to reach | reach before ig + 1) = = o(e e 

k=l 

On the other hand, for mg < ig 

K mg (Tg) < Ei(time to reach ig + 1 | reach ig + 1 before 0) 
+ Ej e+ i(time to reach ig) 

ig oo 1 

i=l k=ig + l 



and for > ig 



oo 



^m g (Tg) = E me (time to reach ig) = V] n Av\m 9 _ Q ( e e/8y 

k=lg + l 



□ 



Proof of Theorem |4.1| (ii). Put 

Sg{5) = inf{t >0:X (t)> 5/9}. 
In view of the previous analysis it is enough to show that for some fixed 

P(m,m)0(9 > S e (5)\T = oo) 1. 
We verify this by showing that for some fixed a G (0, p G (0, j^) 

^{m,m)(re > Sg{ P 9 l - a )\T Q = oo) -)■ 1, (19) 

and 

P m (r e > 5«,(<5)|r«, > ^(p* 1 " )) -> 1. (20) 
14 



According to (11) we have 



\m,m)(Uk = V k ,k = 0,...,6 a |r = oo) 



Note that 



F m (V e - a > P e- a \V k ^0;l<k< 



1, 



since under the condition of non-extinction V n is just a simple random walk 
restricted to the set of positive integers, having a drift that is bounded from 
below by > 0. Combining the last two relations we arrive at (19). 
Finally, ( 20 ) follows from the fact that the probability 

Pj (reach n before 0) = IIj/n n 

with % = p9~ a and n = 5/9 tends to one as 9 — > 0. 

□ 



6 Proof of Theorem 4.2 (ii) 



In the subcritical case s = X/fx lies in (0, 1). To establish the approximation 



formula ( 13 ) we refer to ([8]) which gives 

TO— 1 oo 



E -(^) = E E 



7T 



TO— 1 oo 



J 



j=0 k=j+l 



XklTk 



E 



i=0 k=j+l 



where r,- 



^q^g, and on the other hand, 



TO— 1 oo 7 



IE m (T ) = A- 1 ^ £ 



j=0 k=j+l 



Thus in view of ( 10 ) we have to verify that 



a/e 

E E 

J=0 k=j+l 



22, s k-j 



'i \ v A . /i-A + a , 01 , 

^-^...r^O^—ln-^-^-. (21) 



To prove (21 ) choose arbitrary but fixed small e and large T and consider 



aid j+T 

E E -^-■( 1 -rj---r k -i) 

j=e/9 k=j+l 



a/6 j+T fc ■ 

s r-, v(kO)-vue) 

E E ^•( 1 - e 6 ) + o ^' 



j=e/e k=j+l 
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where a counterpart of Lemma 4.3 was used with 



V(x) = (x + u) In — — x, 

V(x) - V(y) = {x-y) In ^ + {y + /i In — ^ 

/x V y + /j, y + /j, 



v(ke)-v(je) 

Due to the last equality we can replace e » 



with {j^) k 3 and get 



a/0 i+T k _j a/e T i 

Since 



^ s »(l-(l+^/ M )- , ) = r 



J'0 



Z=l 



s I- s{\+ jO/fi)- 1 1 - s n(l - s) +j6' 



to derive (21) it remains to observe that 



a/8 t , 



a/d 



EE ■('-('Wl^Es 



j=e/6 1=1 

where 
and 



j=e/e 



+ Pt(0), 



limsuplimsup |/>r(#)| = 0, 



a/9 



n(l-s)+30 



dx 



In 



{j, — A + a 



— s) + x \i — A + e 



This finishes the proof of (13). 

Next we prove the weak convergence stated in the subcritical case. Fix 
some < a < \. Following the approach of [3], we establish (14) after 
splitting the extinction time tq in two parts 

re = Te t \ + rgp, 



where rg t i is the time for Xg(-) to reach the level a and tq^ is the time 
for the process Xq(-) starting from 9~ a to get absorbed at 0. 
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If Xg(0) = mg and 9mg — > a > 0, then according to [9] the scaled process 
6Xg(-) converges in probability, uniformly on compact time intervals, to the 
deterministic motion x(-) governed by the differential equation 

x'(t) = (X-fi)x(t) -x 2 (t), x(0) = a. (22) 

This equation has an explicit solution 

' n +JL-) e ^--L-. (23) 



x(t) \a // — Ay /i — A 

Solving formally x{t) = 6 l ~ a for the time t required for the deterministic 
motion to reach the low level we find 

(1 - a) h «-■-!■■(<.- + &■ -A)-') 
\i — A 

in probability. Combining ([I]) with ( |11[ ) entails 



which together with (24) give (14). 

The full justification of (24) can be achieved using the approach devel- 
oped in [2] and [3]. It is based on an appropriate integral of the equation 
(22), which in our case is 

, lnfyu — A + z) — lna; + lna — ln(u — A + a) 

h{z,t) = t . (25) 

/U — A 



If x(t) satisfies (23), then h(x(t),t) = and furthermore, x(t — h(z,t)) = z. 
It follows, 

\z-x(t)\ < (/j- X + a)x(t)\h(z,t)\. (26) 



For the rest of the proof we replace a by 9mg in relations (23) and 



( |25| ) defining x(t) and h(x,t). Let z/ e denote the minimal t > such that 
|0JT fl (t) - x(t)| > e, and put He(i) = \h(6X g (t),t)\ so that iZe(O) = 0. 
According to [3] a modified Corollary 1 of Lemma 5 in [2 J gives 

F me (H e (tAv e )>y)<2exp{-ky + tC e (6,k,t)} 

for all positive t and k, where the function C e (0,k,t) can be chosen such 
that for some positive constants Ci, C2, C3 

G (6*, fc, t) = -—— + — — t exp 



(x(t)-e) 2 (x(t)-e) 2 ^\x(t)- 
if we assume that x(t) > e. If furthermore, x(t) — e > C4# 1_a , then 

c(0,m) < c 5 ke 2a - 1 + c 6 k 2 e 2a - 1 e c ' kea . 
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7 Proof of Theorem 4.2 (iii) 



According to ([8j) and (15) we have in the critical case 
eE mg (T d )~ I ^—1. I ^ e B dxdy 



o W Jo V/^+x 



m Jo Vv + X 

— 2/Lt' 



where W(x) = x - (pL + x)ln^. Notice that W(x) = -#(1 + 2r(x)), 



where r(x) — > as x — > 0. It follows, that for any T > 

1 f Ty ^ B ] JJTTy f y 1 w(y)-w{ X ) 

e 9 dxdy 



i^o Jo m Jo Vv + x 



*2 .2 



Jo ^ll + tyfj£ 



T 



— > \i 1 / z 1 / e 2 Z 
Jo Jo 

On the other hand, since for < x < y 

/i + x (y - x) 2 



W{y) - W(x) <(x-y) In 



/i 2(/i + y)' 



we have with C = 



Tv^e y 7o Vf 1 + X 

^ C / 2/ / e 2(M+a)e cte + / e 9 f dx \dy 



T^JIB \J0 Jy/2 
C f°° v 2 f a , f y *-v i n M+y/2 

< — / e s(M+a)e dy + C / y 1 / e 8 m cfedy 



T^P JTy/JS Jy/2 

< — — / e 8 (m+-) dz + Ce 1 



where the last integral is estimated from above by a constant plus 
[ 2 » dy < f 1 < 2y/n + 1 ^ f 4/x 



Tv ^yln(l + ^) ~ J Ty/mr)z ( <z -*) ~ tV6 2 \T^j£ 
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Using a table integral 



oo 

1 



t~—z 







we conclude that 

1 f a JU + V f V 1 W(y)-W(x) 1 /7T\3/2 

— / ; — r u d:nl(i ( - ) £ < r . 



lim sup 

0->O 



f \/fJ> + y f V 1 W(y)-W(x) 1 /7T\ 

/ / , e « array — 

Jo W Jo + x x/l l ^ 2J 



where et — > as T — > oo. 
It remains to observe that 

^t±y.e w ^ $ dy r —^e- w ^/ 9 dx = o(V9). 

m Jo \Jv + x 

Remark. Our approximations for the mean extinction time are specific to 
the population model we study. These should be compared with similar 
calculations performed in a more general setting by [3], where, however, 
strict justifications of some important steps are missing. 
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